First, we’ll make the main sampling map.
mapbox <- c(-76.14, 43.14, -75.7, 43.3) # specify map boundaries
oneida_lake <- get_map(location = mapbox, source = "stamen",
maptype = "terrain", zoom = 13)
# Could pull up google maps if you have a license
cols <- c(brewer.pal(8, "Dark2"),"#386CB0")
oneida_lake_map <- ggmap(oneida_lake) +
geom_point(data=all_sampling_points,
mapping=aes(x=long, y=lat, color=gear, shape=gear),
size=3, stroke = 1.2) +
scale_color_manual(values=cols[c(4,9,6,3,1)]) +
scale_shape_manual(values=c(17,5,3,16,6)) +
geom_line(data=ef_dat, aes(x=long, y=lat, group=Description),
color="white", size=0.5) +
theme(axis.text.x = element_blank(),
axis.text.y = element_blank(), axis.ticks = element_blank(),
rect = element_blank(), axis.title.y=element_blank(),
axis.title.x=element_blank()) +
ggsn::scalebar(location = "bottomleft", x.min = -76.13, x.max = -75.7,
y.min = 43.145, y.max = 43.25,
dist = 5, dist_unit="km", transform = TRUE,
model = "WGS84", height = 0.05,
st.dist = 0.05, st.bottom=FALSE)